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Atmospheric  propagation  of  partiaily  coherent  radiation* 

J.  Carl  Leader 

McDonnell  Douglas  Research  Laboratories,  McDonnell  Douglas  Corporation,  St.  Louis,  Missouri  63166 

(Received  17  August  1977) 

An  extended,  Rayleigh-Sommerfeld  integral  method  is  used  to  derive  expressions  for  the  mutual  co¬ 
herence  function  and  radiation  intensity  derived  from  a  planar,  partially  coherent  source  propagating 
through  the  atmosphere.  The  derived  results  reduce  to  previous  results  for  (i)  coherent  radiation  propa¬ 
gation  in  the  atmosphere  and  (ii)  the  relations  relating  the  far-field  intensity  angular  distribution  and 
the  source  coherence  for  a  partially  coherent  source  in  vacuo,  '  A  mathematical  description  of  the  pre¬ 
dicted  results  in  terms  of  the  vacuum  distribution  and  scattering  functions  (related  to  the  Fourier-trans¬ 
formed  two-source  mutual  coherence  function)  is  permitted  by  this  development.  Analytical  results  are 
calculated  for  a  homogeneous  atmosphere  and  a  source  coherence  that  simulates  a  laser- illuminated 
rough  surface.  The  effective  far-field  range  is  determined  by  the  source  size,  wavelength,  and  source  co¬ 
herence  length.  The  phase  of  the  calculated  mutual  coherence  function  is  determined  by  the  field-point 
separation  for  off-axial  propagation  directions.  Numerical  results  for  the  amplitude  and  phase  coher¬ 
ence  lengths  are  calculated  and  illustrated  as  a  function  of  the  source  size,  source  coherence  length,  prop¬ 
agation  angle,  range,  and  refractive -index  structure  constant. 


INTRODUCTION 

Numerous  investigators  have  studied  the  coherence  of 
radiation  propagating  in  the  turbulent  atmosphere  using  a 
variety  of  approaches  including  difference  equations,^ 
transport  methods, ^  the  Markov  approximation,^  ladder  ap¬ 


proximations  of  the  Bethe-Saltpeter  equation,^  and  the  Hu- 
ygens-Fresnel  principle.*^  Reran®  has  provided  an  excellent 
account  of  the  consensus  of  these  various  approaches;  All  of 
the  aforementioned  studies,  however,  have  been  principally 
concerned  with  the  propagation  of  initially  coherent  radiation 
because  of  its  relevance  to  laser  propeigation  phenomena.  Lee 
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etalJ  have  examined  the  propagation  of  perfectly  incoherent 
radiation  in  a  turbulent  medium,  using  the  extended  Hu- 
ygens-Fresnel  principle,  to  determinie  atmospheric  effects  on 
speckle  phenomena.  None  of  the  previous  studies  have 
treated  in  detail  the  problem  of  atmospheric  propagation  of 
radiation  derived  from  a  partially  (spatially)  coherent  source. 
The  purpose  of  this  paper  is  to  examine  the  general  case  of 
partially  coherent  radiation  propagation  in  a  random  rnedi- 
um. 

The  approach  used  in  this  study  is  essentially  an  extended 
Rayleigh-Sommerfeld  integral  method.  This  technique  is 
analagous  to  the  extended  Huygens-Fresnel  method®;  how¬ 
ever,  it  facilitates  the  analysis  of  the  behavior  of  radiation 
coherence  as  a  function  of  the  propagation  angle.  Techniques 
previously  employed  to  study  coherent  radiation  propagation 
use  coordinates  appropriate  for  (approximate)  rectilinear 
propagation.  However,  partially  coherent  radiation  propa¬ 
gates  throughout  an  angular  spectrum.®  Therefore,  polar 
coordinates  are  more  appropriate  than  Cartesian  coordinates 
to  study  partially  coherent  radiation  propagation.  Impor¬ 
tantly,  the  derived  results  relax  to  the  Huygens-Fresnel  result 
for  the  special  case  of  a  coherent  source,  demonstrating 
equivalence  with  the  previous  coherent  source  analyses,^® 
Significantly,  they  also  relax  to  the  far-field  coherence-angular 
distribution  relations®* characterizing  a  partially  coherent 
source  for  the  special  case  of  vacuum  propagation.  A  simple 
physical  picture  of  the  atmospheric  effects  in  terms  of  scat¬ 
tering  functions  is  an  additional  benefit  provided  by  this  ap¬ 
proach. 

To  illustrate  the  effects  of  range,  propagation  angle,  source 
size  and  coherence,  and  strength  of  turbulence  on  the  spatial 
coherence  (second  moment)  predicted  by  this  analysis,  sim¬ 
plifying  approximations  are  made  which  facilitate  closed-form 
solutions.  The  phase  of  the  predicted  mutual  coherence 
function  (MCF)  changes  with  field-point  separation  at  off- 
normal  propagation  angles.  This  effect  has  not  been  pre¬ 
viously  predicted  because  of  assumed  axial  propagation  (in¬ 
cluding  the  speckle  analysis  of  Lee  et  al.  '^).  Numerical  results 
for  the  MCF  amplitude  and  phase  coherence  lengths  are  cal¬ 
culated  and  depicted. 

THEORETICAL  FOUNDATION 

A  monochromatic  (spatially)  partially  coherent  source  is 
assumed  that  occupies  a  portion  of  the  x~y  plane  (Fig  1). 
Physical  embodiments  of  such  a  source  are  laser-illuminated 
diffusers  or  rough  surfaces.  Restricting  the  analysis  to  the 
scalar  problem^^  for  simplicity  and  ease  of  comparison  with 
previous  work  and  denoting  the  vectors  r  =  (x,  y,  0)  and  P  = 
U,  y,  z),  and  the  wave  number  ky  the  scalar  field  at  P  (resulting 
from  the  source  propagating  through  a  turbulent  medium)  is 
given  by 

^  J  J_- 

X  exp[i/2|  r  -  P|  +  \p{Py  r)],  (1) 

where  xj/  (P,  r)  represents  the  phase  and  log-amplitude  vari¬ 
ation  from  the  vacuum  solution  for  a  spherical  wave  propa¬ 
gating  from  r  to  P.  Equation  (1)  follows  from  an  extension 
of  the  ‘‘self-consistent”  Rayleigh-Sommerfeld  solution^®*^^ 
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for  the  scalar  field  at  P  resulting  from  the  aperture  field  U (r). 
The  extension  of  this  solution  to  turbulent  medium  propa¬ 
gation  consists  essentially  of  the  additional  complex  phase 
term  ^(P,  r)  following  identical  arguments  presented  in  Refs. 
8  and  15.  Assuming  that  the  source-field  distribution  is  rel¬ 
atively  large  with  respect  to  a  wavelength  and  restricting  the 
analysis  to  propagation  distances,  R-  (x2-|-y2  +  22)1/2 
are  much  greater  than  a  typical  source  dimension,  Eq.  (1)  can 
be  written 


[/t>xb(p)s_£2!^  r  r  “  dWWc'M r-P|+^(r.P) 
27rjR  ♦/— 00 

(2) 

to  an  excellent  approximation,  where  0  is  the  angle  subtended 
between  P  and  the  z  axis.  Denoting  the  average  over  the 
ensemble  of  all  source  and  turbulent  medium  configurations 
by  angle  brackets  (  > ,  the  field-point  MCF  can  be  written 

rturb(p^p/)  =  (f/turb(p)(7*turb(p/)) 


cos6  cosd'k^ 


<xt  ^  ^00 

X  F(r,  r',  P,  PO  exp[ikL(ry  r\  P,  P%  (3) 


where 

r,(r,  r')  =  {U(r)UHr'))y  (4) 

F(r,  r',  P,  F)  =  <exp[^(r,  P)  +  ^*(r',  P')]>,  (5) 

L(r,r',P,PO  =  |r-P|-|r'-F|,  (6) 

and  the  order  of  integration  and  ensemble  averaging  has  been 
interchanged.  Note  that  the  averages  of  the  source-field  and 
medium  propagation  variables  are  separable  because  these 
are  presumably  independent  random  processes.  Denoting 
unit  vectors  s  and  s'  in  the  respective  directions  P  and  P',  the 
following  sum  and  difference  coordinates  are  defined: 
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r+  =  (r  +  rO/2  and  r“  =  r  ~  r',  (7) 

P+  =  (P  +  P0/2  and  P-  =  P-P',  (8) 

=  (s  +  sO/2  and  ,s~  =  s—  s\  (9) 

R+={R  +  R')/2  and  R- =  R  -  R\  (10) 

6+ =  {6  +  d')/2  and  d  =  6-  0\  (11) 


Appendix  A  demonstrates  that  the  phase  factor  L(r,  r',  P,  P') 
is  approximately  given  by 

L(r,  r',  P,  PO  s  (r  •r+/i2+  -  r”  -  r'^ -5“  +  R~),  (12) 

The  ensemble-averaged  medium  propagation  factor  F(r,  r', 
P,  P')  is  recognized  as  the  two  (point) -source  spherical- wave 
mutual  coherence  function.  Assuming  isoplanicity,  the 
variable  dependence  of  the  two-source,  spherical  wave  MCF 
is 

F(r,r^P,P')-F(r~  P-;/?+).  (13) 

The  field  point  MCF  therefore  can  be  written 

inurb(p+  p~)^  _ r _ _ 

(2irfl  +  )2  [1  -  (R-/2R+)^ 

X  j  c(2r+d2r-  r,(r+  r-)F(r-  P";  F+) 

X  exp[t/e(r“  •  r+//?+  “  r“  •  -  r+  •  s'")],  (14) 

assuming  small  field-point-difference  angles  5.  Defining  the 
function 


G(\,  P-)  =  J"”  d2(Ar-)e-'(*'-)-«G(r-  P"),  (21) 


and 


hi,  P")  =  J"  c(2(/jr-)e-‘(*>-)-fF(r-  P").  (22) 

Because  the  vector  f  corresponds  to  a  unit  vector  direction, 
the  Fourier-transformed  generalized  two-field -point  source 
coherence  can  be  expressed 

r,-,  T.-1  I'  -  9-m'n 

’  it  (i-il-i?) 

X  P")  (23) 

using  Eq.  (20)  for  the  vacuum  MCF.  Thus  the  turbulent 
medium  MCF  can  be  expressed  in  terms  of  the  vacuum  MCF 
as 


rturb(p+,  p-) 


COS^^'*' 

7^ 


xx: 


d2| 


rvac(F|,  P- 

(1  -  -  x5) 


) 

-F(s+ 


P") 


(24) 


using  Eqs.  (17),  (22),  (23)  and  the  convolution  theorem  of 
Fourier  transformation.  Note  that  for  vacuum  propaga¬ 
tion 

^(t  P-)  I  vacuum  =  (27r)25(|),  (25) 


using  Eq.  (19),  so  that 


Gir-,  P-)  . 

=  IM  J'  J'“  d2r+r,,(r+  (15) 

as  the  generalized,  two-field-point,  source  coherence  (the 
rationale  for  this  appelation  is  provided  below),  where 


A=J^d2r+  (16) 

is  the  source  area;  the  field-point  MCF  is 

rturb(p+  P-)  =  eikR-  cos2e+ 

(2ir)2[l-(F-/2/?+)2] 

X  J*  J^”  d2(Ar-)e-‘(*‘-)-‘«''G(r-  P-)F(r-  P-),  (17) 

where 

fi3  =  A(F+)~2  (18) 


is  the  solid  angle  subtended  by  the  source  at  range  R'^,  It  is 
noted  that  for  vacuum  propagation 

F(r-P-)  I  vacuum  =  1  (19) 


SO  that  the  field-point  MCF  reduces  to 
\  cos^^'*' 

pvac(p+  p— ^  ®  '' 


(27r)2[l-(F-/2F+)]2 

^  S  S-  P-)  (20) 


for  vacuum  propagation.  Spatial  source-plane  Fourier 
transforms  of  the  generalized  two-fleld-point  source  coherence 
and  the  two-source  spherical-wave  MCF  are  now  defined. 


r‘“'-b(P+p-)  I  vacuum  =  C0s2r 

^  /  XI  =  r™^(P+,  P-),  (26) 

providing  a  check  on  the  self-consistency  of  Eq.  (24). 

Partial  validation  of  the  above  equations  is  obtained  by 
examining  relevant  limiting  cases.  Clearly,  as  the  source 
becomes  perfectly  coherent,  the  radiation  propagates  pri¬ 
marily  along  the  z  axis.®  Thus,  field-point  vectors  (p  and  pO 
that  are  transverse  to  the  z  axis  can  be  defined  to  parameterize 
the  field  point  separations.  The  unit  direction  vectors  s  and 
s'  are  then  given  by 

s  ^  p/2  +  k 

and  ,  (27) 

5' s  pV2'  +  k, 

where  ^  is  a  unit  vector  along  the  2  axis.  Substitution  of  Eq. 
(27)  in  Eq.  (14)  then  yields^^ 

n^’-Kp,  pO  = /22/(27r2)2 

X  r  f'”  f  f“d2r+d2r-i/(r+-hr-/2) 

X  f7*(r+  -  r“/2)F(r~,  p") 

X  exp[ik{r+  -  p+)  •  (r"  -  p”)],  (28) 

where 


P+  =  (p  +  p')/2,  p-  =  p-p', 

and  the  coherent  form  for  the  source  MCF^'^  has  been  used 
together  with  the  fact  that  cos^'*’  ^  1  and  the  assumption  that 
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z-z'.  Equation  (28)  is  recognized  as  the  turbulent  medium 
MCF  predicted  by  the  extended  Huygens-Fresnel  principle 
for  a  coherent  source/'’ 


Examining  the  far-field  propagation  of  the  radiation  in¬ 
tensity  in  vacuo,  Eq.  (14)  reduces  to 

^''“(P+)|far-field  =  ^t-'’(P^O)|  vacuum, far-field  ~ 

X  J'  J"”  J  J*”  d^r^d^r-  r,(r+,  (29) 

Assuming  that  the  source  is  statistically  homogeneous  (in 
the  same  sense  described  in  Ref.  9),  Eq.  (29)  reduces  to 

jvac(p+)  I  =  A  ts  (ksl ),  (30) 

where 


f*(^sl)  =  f_fl  d^r- e-‘*'-^*r,(r-),  (31) 

the  Fourier  transform  of  the  source  mutual  coherence,  where 
Sx  is  the  projection  of  s  on  the  x-y  plane.  Equation  (30)  is 
equivalent  to  the  relationship  derived  by  Wolf  and  Carter^ 
relating  the  far-field  intensity  angular  distribution  and  the 
source  coherence.  Furthermore,  an  equivalent  statement  of 
the  far-field  van-Cittert-Zernike  theorem  can  be  derived  using 
Eqs.  (14)  and  (30),  employing  the  same  procedures  described 
in  the  reference  cited  in  Ref.  12. 


Returning  to  the  problem  of  turbulent  medium  propaga¬ 
tion,  it  is  noted  that  the  radiation  intensity  in  the  turbulent 
medium  is  given  by  [Eq.  (24)], 


Jturb(p+)  =  nurb(p+^  0)  == 


cos^^'*' 

(27r)2 


X 


Jvac(j^+|) 

a-kl-  ey) 


A2f(l+  - 1). 


(32) 


where 


jGr(^)  =  J'  J'”  c/2(ftr-)e--<**-')-ff(r-  0)  (33) 

is  the  Fourier-transformed  mutual  coherence  function  for  a 
spherical  wave  propagating  in  the  turbulent  medium. 
Equation  (32)  lends  itself  to  a  simple  physical  interpretation 
since  it  is  exactly  the  same  form  that  results  from  postulating 
a  scattering  function  that  describes  the  angular  deviation  of 
the  vacuum  radiant  intensity  resulting  from  propagation  in 
the  turbulent  medium.  The  total  intensity  in  any  given  di¬ 
rection  is  then  the  convolution  of  the  vacuum  intensity  dis¬ 
tribution  with  the  postulated  scattering  function.  Clearly, 
in  this  case,  the  scattering  function  is  given  by  the  Fourier- 
transformed,  spherical-wave  MCF  divided  by  the  factor  (1  - 
Analogously,  Eq.  (24)  expresses  the  turbulent  me¬ 
dium  MCF  in  terms  of  the  vacuum  MCF  angular  distribution 
(for  a  fixed  field-point  separation  vector  P“)  convolved  with 
a  pseudo-scattering  function  that  describes  the  average  vector 
angular  deviation  of  the  vacuum  MCF  (again  with  a  fixed 
field-point  separation  vector  P")  resulting  from  propagation 
in  the  turbulent  medium.  In  this  case  the  pseudo-scattering 
function  is  given  by  the  Fourier-transformed,  two-source, 
spherical  wave  MCF  divided  by  the  factor  (1  -  To 

complete  the  physical  interpretation  of  the  above  equations, 


it  is  noted  that  the  far-field  expression  for  G(r  ,  P  )  [Eq.  (15)] 
with  zero  field -point  separation  is 

G(r-  0)|far-ficid  =  t  J  d2r+r,(r+  r"),  (34) 

which  is  the  spatial  average  of  the  source  mutual  coherence 
function.  Thus,  the  general  function  G(r“,  P”)  is  the  gen¬ 
eralization  of  the  average  source  mutual  coherence  that  results 
from  considering  separate  field  points  and  near-field  ranges. 
Hence  the  appellation  of  the  generalized,  two-field-point 
source  coherence  for  the  function  G  (r~,  P“)  results.  Equation 
(17)  provides  a  mathematical  statement  of  the  propagation 
of  the  (generalized)  source  coherence  in  a  turbulent  medium. 
The  Fourier-transform  relations  that  have  been  noted  permit 
an  interpretation  of  this  expression  in  terms  of  scattering 
functions. 

CALCULATIONS 

To  perform  specific  calculations  of  the  MCF  obtained  from 
a  partially  coherent  source  propagating  through  the  turbulant 
atmosphere,  it  is  necessary  to  specify  the  generalized,  two- 
source  spherical- wave  MCF,  F(r“,  P“),  or  equivalently,  the 
pseudo-scattering  function  J^(|,  P  ” ) .  If  the  statistics  of  the 
complex  phase  \l/  are  Gaussian,  then  the  two-source  spheri¬ 
cal-wave  MCF  is  given  by^® 

F(r-  P-)  =  exp  [-^D^(r-  P-)  j,  (35) 

where 

D^(r-  P-)  =  D^(r-  P")  +  Ds(r-  P"),  (36) 

the  sum  of  the  log-amplitude  (D^)  and  phase  (0$),  two-source 
spherical-wave  structure  functions.  Assuming  a  Kolmogorov 
turbulence  spectrum  and  weak  turbulence  the  complex  phase, 
two-source,  spherical-wave  structure  function  is  given 

byl9.20 

DJr-P)  =  2,91^2/2+  C\tCUR-^t)\t^P(P-) 

%/0 

+  (l-t)^^(r-)|5/3,  (37) 

where  C^(R)  is  the  refractive  index  structure  constant,  radial 
changes  in  the  complex  phase  \p  are  ignored  (these  are  small 
for  small  R"),  and  the  operator  P(x)  projects  the  vector  x  onto 
a  plane  that  is  transverse  to  the  P"*"  direction.  Equation  (37) 
assumes  that  the  absolute  values  of  the  projected  vector  dif¬ 
ferences  |"P(r“)|  and  lP{P~)\  are  much  less  than  a  Fresnel 
zone  length,  (27rR'^ Although  the  integral  of  Eq.  (37) 
cannot  be  integrated  exactly,  approximate  expansions  can  be 
obtained  and  used  with  Eq.  (35)  and  the  preceding  develop¬ 
ment  to  perform  numerical  calculations  of  the  MCF  of  a 
partially  coherent  source.  However,  these  procedures  are 
complicated  and  permit  little  physical  understanding  of  the 
propagation  characteristics  as  a  function  of  the  variables  of 
interest.  Analytical  calculations  of  the  MCF  are  permitted 
by  approximating  the  five-thirds  power  of  the  projected  vector 
sum  in  Eq.  (37)  by  the  square.  Normalizing  the  resultant 
two-source  spherical-wave  MCF  to  the  correct  five-thirds  law 
spherical -wave  coherence  length  at  the  e~^  point  yields 
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F(r-  P-)  s  exp  Spo"^  dt\tP{P-) 

+  (l-t)P(r-)pj  ,  (38) 

SO  that 

M(p)  ^  F[0,  p  =  |P(P-)|]  -  (39) 

where 

Po  =  [%  X  L455/2  2i?+C2]-3/5  (40) 

assuming  homogeneous  turbulence.  Appendix  B  provides 
a  calculation  of  the  MCF  of  a  uniform  coherent  source  using 
the  approximations  of  Eqs.  (38)-(40)  and  the  preceding  de¬ 
velopment.  It  is  noteworthy  that  the  plane-wave  coherence 
length  predicted  by  this  approximation  (0.58po)  is  in  reason¬ 
able  agreement  with  the  five- thirds  law  prediction  (0.56po). 
Calculations  are  also  provided  in  Appendix  B  for  the  Fourier 
transform  of  the  uniform  source  mutual  coherence  in  terms 
of  a  dimensionless  spatial  frequency  spectrum  V  and  a  di¬ 
mensionless  spatial  frequency  x*  The  predicted  spectrum  of 
the  mutual  coherence,  resulting  from  the  quadratic  approxi¬ 
mation,  Eqs.  (38)-(40),  is  plotted  in  Fig.  2  along  with  experi¬ 
mental  data  of  Artem’ev  and  Gurvich^^  and  three  different 
calculations  of  a  modified  Kolmogorov  spectrum^i  corre¬ 
sponding  to  three  different  values  of  the  inner  scale  of  tur¬ 
bulence.  Significantly,  the  approximate  quadratic  spectrum 
is  within  1  dB  of  the  experimental  data  over  the  range  of 
frequencies  plotted  and  is  closer  to  the  measured  data  than 
two  of  the  modified  Kolmogorov  spectra  for  the  highest 
measured  spatial  frequencies.  Because  the  quadratic  ap¬ 
proximation  is  in  reasonable  agreement  with  available  MCF 
spectral  data,  it  is  employed  for  subsequent  calculations  to 
obtain  closed-form  analytical  expressions  for  the  MCF  of  ra¬ 
diation  propagating  from  partially  coherent  sources.  It  is 
expected  that  although  detailed  predicted  behavior  may  de¬ 
viate  somewhat  from  that  obtained  from  the  five -thirds  law 
MCF,  the  predicted  trends  with  respect  to  the  various  vari¬ 
ables  of  interest  should  be  the  same  as  those  obtained  from 


FIG.  2.  Experimental  spatial  spectra  of  the  coherence  function  obtained 
by  Artem’ev  and  Gurvich^”'  (open  circles)  compared  with  predictions  derived 
from  a  modified  Kolmogorov  turbulence  spectrum  (solid  lines)  and  the 
quadratic  approximation  (heavy  solid  line)  discussed  in  Appendix  B. 
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the  more  accurate  description  of  the  atmospheric,  spherical- 
wave  structure  function. 

The  form  of  the  source  mutual  coherence  function  that  will 
be  used  for  the  subsequent  investigation  is 

r,(r+  r-)  =  /,7rexph(|r-|2/p2+  jr+p/L^)].  (41) 

Equation  (41)  is  particularly  suitable  for  this  investigation 
because  when 

Pg  ~  2Zy, 

r,(r,  ro  =  /.TT  exph(2L)-2(|r|2  -f  |r'|2)],  (42) 

which  is  the  form  of  the  source  coherence  for  a  coherent 
Gaussian  beam.  However,  when  <  2L,  the  source  coher¬ 
ence  is  similar  to  that  obtained  from  specular  glints  from  a 
'  laser-illuminated  rough  surface  with  an  equivalent  rms  surface 
slope  given  by 

(43) 

Thus  variation  of  the  paraixieter  pe  simulates  laser  reflectance 
(TEMoo)  from  a  surface  with  a  variable  surface  roughness. 
The  similarity  of  the  source  MCF  with  that  of  the  rough  sur¬ 
face  is  demonstrated  by  noting  that  the  far -field  vacuum  in¬ 
tensity  distribution  associated  with  the  assumed  source  co¬ 
herence  is  [Eq.  (29)] 

J""^(^)|far-fielcl  =  COs2P+(/zp,)2 

X  exp[— (/2pe/2)2  sin2^]  (44) 

using 

A  =  y  J'”  d^r+  =  wL^.  (45) 

Recalling  the  Fourier-transform  relationship  between  the 
far-field  intensity  angular  distribution  and  the  source  co¬ 
herence  [Eq.  (30)  and  Refs.  9  and  11]  and  noting  that  the  an¬ 
gular  distribution  of  the  radiation  intensity  scattered  by 
specular  glints  of  a  rough  surface  is22 

«/(^'^)  I  specular  glints  ~  (secWs2)  exp(-tan2p/s2),  (46) 

where  s  is  the  rms  surface  slope,  the  stated  equivalence  [Eq. 
(43)]  is  established  in  the  small-angle  approximation. 

Before  proceeding  with  calculations  of  the  turbulent  me¬ 
dium  propagation  of  radiation  derived  from  the  assumed 
partially  coherent  source  [Eq.  (41)],  it  is  instructive  to  examine 
vacuum  propagation  of  the  mutual  coherence.  The  vacuum 
MCF  predicted  by  Eq.  (20)  for  the  source  MCF  of  Eq.  (41) 
is 

X  exp  (^)  (R+)-i(/3sJ  -  cos0+)] .  (47) 

where 

/3  =  o)  sin0+,  (48) 

Sx  =  sin^'*'  sin</>,  (49a) 

Sy  =  sin^'*'  cos(^,  (49b) 
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mRfr/R^)  =  [1  +  (Rf,/R^)T\  (50) 

Rfi  =  kLpe/2,  (51) 

to  is  the  azimuthal  difference  angle  between  the  two  field-point 
coordinates,  and  (j)  is  the  azimuthal  angle  subtended  between 
the  y  axis  and  the  projection  of  the  vector  on  the  x-y  plane. 
Note  that  the  radiation  intensity  predicted  by  Eq.  (47)  (i.e., 
d  =  ^  =  R~  =  0)  approaches  the  predicted  far-field  intensity 
when  /?+»/?  ff.  Therefore  R ff  is  interpreted  as  the  equivalent 
far-field  range  for  a  partially  coherent  source  with  a  Gaussian 
coherence  function.  This  interpretation  is  supported  by  the 
following  observations.  When  the  source  coherence  length 
is  equal  to  the  source  dimension  (2L),  the  equivalent  far-field 
range  is  (R((  =  2TrL^/X),  within  a  factor  of  1.27  of  the  far-field 
range  customarily  predicted  for  a  uniformly  illuminated  ap¬ 
erture  [2{2L)y\],  However,  when  the  coherence  length  of  the  ‘ 
source  is  on  the  order  of  the  wavelength,  the  equivalent  far- 
field  range  decreases  to  the  order  of  the  source  dimension 
(ttL).  It  is  well  known  that  valid  measurements  of  the  far- 
field  angular  distributions  of  incoherent  sources  can  be  per¬ 
formed  at  close  range.  The  equivalent  far-field  range  pa¬ 
rameter  Rff  critically  affects  the  range  dependence  of  the 
vacuum  MCF  and  the  turbulent  medium  MCF  calculated 
below.  As  the  far-field  range  is  exceeded,  the  MCF  predicted 
by  Eq.  (47)  approaches  the  prediction  of  the  far-field  van- 
Cittert-Zernike  theorem  for  the  assumed  source  irradiance 


distribution.  Because  the  phase  difference  of  the  field  points 
[predicted  by  Eq.  (47)]  increases  with  increasing  source  co¬ 
herence  length,  the  phase  portion  of  the  predicted  MCF  re¬ 
sults  from  deterministic  beam-spreading  effects. 

Appendix  C  provides  a  derivation  of  the  generalized  two- 
source  spherical-wave  mutual  coherence  function  F(r“,  P“, 
/?■'■)  [using  the  approximations  of  Eqs.  (38)-(40)]  and  the 
pseudo-scattering  function  F(l,  P“;  R*^).  The  derived  re¬ 
sults  are 

F(r-,  P-)  =  exp  I-  {R+/ponib  + 

-l-■(c  -  (5)2-hc6  -  5/5]!,  (52) 


b  =  x-/R+,  c=y- cose+/R+,  (53) 

Fil  P-  F+)  =  7r(/2po)2  sec0+  exp-[%(F+/po)2(/52  +  52)] 

X  exp  y4{kpo)y^x  +  sec26>]} 

.X  exp[i{kR'^/2)(^xl3  -  sec(9)].  (54) 

Employing  Eq.  (52)  with  Eqs.  (17)  and  (15)  [or  equivalently 
Eq.  (54)  with  Eqs.  (24)  and  (47)]  yields  the  following  MCF 
resulting  from  the  assumed  partially  coherent  source  [Eq.  (41)] 
propagating  through  a  (homogeneous)  turbulent  atmo¬ 
sphere: 


pturb(p-l-^  p-)  = 


Jvac(^-h)|^^r.field  /  oo  \ /R^\H .  ^  n-^~\  I 

Xexp(-62|^— +  (— )%os20+[l-(Q-l)Q'-i]J) 

X  exp  1  (^')"  [sin2^+ -  (sJ)2Q-i  -  (s;)2Q'-i]) 

X  exp j  [  (P)%  1  £)”]9-  - [  (ff)%  1  (^)*] 


Q  =  [l  +  (Fff/F+)2+(p,/po)2] 


Q'  -  [1  +  (Fff//?+)2  -I-  (pe/po)^  COs2^+].  (56b) 

It  can  be  verified  by  inspection  that  Eq.  (55)  reduces  to  the 
vacuum  result  as  the  atmospheric  turbulence  goes  to  zero  (i.e., 
Po  — ^  °° ).  It  is  noteworthy  that  the  ratio  of  the  source  coher¬ 
ence  length  to  the  atmospheric  coherence  length  (p^/po)  is  a 


critical  parameter  in  determining  the  atmospheric  effects  of 
partially  coherent  radiation  propagation.  This  observation 
is  consistent  with  the  analysis  of  Dunphy  and  Kerr ,2*^  who 
noted  the  critical  nature  of  the  ratio  (L/po)  for  the  propagation 
of  coherent  radiation.  To  facilitate  analyses  of  the  effects  of 
various  parameters  on  the  predicted  spatial  coherence,  it  is 
convenient  to  normalize  Eq.  (55)  to  the  field-point  intensities 
(i.e.,  construct  the  degree  of  coherence  function),  choose  an 
azimuthal  plane  of  scattering  (0  =  0),  and  ignore  the  relatively 
uninteresting  azimuthal  dependence  (i.e.,  choose  /5  =  0).  The 
degree  of  coherence  is  then 


^turb(p-f-^  P-)  = 


rturb(p^  PO 


rturb(p+^  P-) 


[pturb(p^  p)pturb(p/^  p^)]l/2  0) 

_l_  cos20+[l  -  (Q  -  l)Q^“i]|^  exp  |-f(5F''')5  sin^  cosO  ~ 


for  small  field-point  separation  angles  <5, 

A  further  simplification  results  from  defining  the  degree 
of  coherence  as 

yid+,  R+-,  d,R-~)  =  T(r,  F+  5) 

X  -  {R-/2R+)^-\  (58) 


where  the  well-known  radial  dependence  of  the  degree  of  co¬ 
herence  has  been  separated  from  the  more  interesting  am¬ 
plitude  and  phase  dependence.  Further,  defining  a  magni¬ 
tude  correlation  angle  de  via  the  equation 

T{d+,  F+  5e)  =  c-iT(0+  F+  0)  =  (59) 
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and  a  {tt/A)  phase  correlation  angle  5^/4  via  a  similar  relation, 
i.e., 

a{6'^,  =  x/4,  (60) 

the  respective  magnitude  and  phase  correlation  angles  follow 
from  Eq.  (57),  i.e.. 


+  ^  cos20+(l  -  (Q  -  l)Q'-i]p^^  (61) 


6W4  =  -(fei?  +  ) 


Equations  (61)  and  (62)  define  the  and  7r/4  points,  re¬ 
spectively,  of  the  amplitude  and  phase  angular  field  correla¬ 
tion.  Corresponding  coherence  lengths  are  defined  by  the 
equations 

A,  =  (63) 


K/a  K/aR'^^  (64) 

Several  observations  can  be  made  from  inspection  of  the 
above  equations.  First,  the  phase  coherence  length  becomes 
infinite  for  the  axial  and  grazing  propagation  directions  (i.e., 
the  MCF  phase  remains  constant  regardless  of  the  field-point 
separation).  This  behavior  is  consistent  with  all  previous 
studies  of  coherent  radiation  propagation  since  phase  effects 
have  not  been  previously  noted.  In  general,  however,  a  fi¬ 
nite-phase  correlation  length  is  predicted  from  this  analysis 
because  of  beam  spreading  effects.  Although  the  predicted 
phase  coherence  length  (at  nonaxial  directions)  increases  as 
the  range  greatly  exceeds  the  far-field  range  and  the  atmo¬ 
spheric  coherence  len^h  greatly  exceeds  the  source  coherence 
length,  these  are  generally  incompatible  simultaneous  re¬ 
quirements  for  atmospheric  propagation  unless  the  source  is 
completely  incoherent  (i.e.,  b  correlated).  Because  perfectly 
incoherent  sources  are  prohibited, finite-phase  coherence 
lengths  are  predicted  at  nonaxial  directions  in  the  atmosphere. 
The  predicted  angular  dependence  of  the  amplitude  coherence 
length  is  also  worthy  of  note.  Clearly,  as  grazing  angle 
propagation  (90°)  is  approached,  the  effects  on  the  source 
coherence  length  resulting  from  the  finite  coherence  length 
are  eliminated  because  of  the  cos^^"*"  factor  multiplying  these 
terms.  The  degree  of  coherence  approaches  that  predicted 
for  a  spherical  wave  [Eq.  (39)]  in  the  polar  plane,  i.e.,  y  ~ 
exp[-(/t!'*’5)V/Oo]j  as  the  grazing  geometry  is  approached.  The 
physical  reason  for  this  phenomenon  is  that  the  projected 
source  area  approaches  a  line  source  as  the  grazing  angle  is 
approached  and  therefore  becomes  perfectly  coherent  in  the 
transverse  direction  (i.e.,  the  polar  plane).  Thus,  a  spherical 
wave  MCF  functional  form  is  expected  in  the  polar  plane  since 
the  source  is  effectively  a  point  source  in  this  dimension. 

To  further  explore  the  variable  dependence  of  the  MCF 
predicted  by  the  above  analysis,  three-dimensional,  computed, 
graphical  results  are  employed.  Calculated  values  of  the 
amplitude  and  phase  coherence  lengths  are  plotted  as  the 
ordinate  variation,  while  various  combinations  of  the  propa- 


J.  Opt.  Soc.  Am.,  Vol.  68,  No.  2,  February  1978 


0  10  20  30  40  50  60  70  80  90 


6  (deg) 


6  (deg) 

FIG.  3.  Calculated  values  of  the  amplitude  (a)  and  phase  (b)  coherence 
lengths  as  functions  of  the  propagation  angle  (0+)  and  range  {R'^).  Fixed 
parameters  are  X  =  10“®  m,  =  10“''®  m”^^^  L=  10"''  m,  and  Pe  = 
10“'^  m. 

gation  parameters  are  used  as  independent  variables  along  the 
X  and  y  axes.  Representative  optical  values  of  the  fixed  in¬ 
dependent  variables  have  been  chosen  for  the  subsequent  il¬ 
lustrations. 

Figure  3  illustrates  the  range  and  angle  dependence  of  the 
amplitude  and  phase  coherence  lengths  [Eqs,  (63)  and  (64)] 
for  a  relatively  incoherent  source  {pe  =  lOX)  and  moderate 
turbulence  (Cl  =  10”i^  m  x  =  iq-s  m).  The  aforemen¬ 
tioned  dependence  of  the  apparent  source  coherence  on  the 
propagation  direction  is  clearly  illustrated  by  the  amplitude 
coherence  length  dependence  on  the  range  and  direction  of 
propagation  [Fig,  3(a)].  It  is  apparent  that  the  amplitude 
coherence  length  is  dominated  by  the  finite  source  coherence 
(i.e.,  speckle  dominated)  for  axial  propagation  [note  that  Ag 
is  approximately  proportional  to  (R'^)%  whereas  atmospheric 
effects  begin  to  interact  and  dominate  as  the  grazing  propa¬ 
gation  angle  is  approached  [thus  Ag  po  iR'^)~^^^].  The 
phase  coherence  length  is  plotted  only  from  5°-85°  because 
it  goes  to  infinity  at  the  endpoints  of  the  angular  distribution 
(0°  and  90° ) .  The  angular  dependence  of  the  phase  coherence 
length  is  relatively  unaffected  by  the  range  for  the  chosen  fixed 
independent  variables. 

To  illustrate  the  effects  of  the  finite  source  coherence  length 
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FIG.  4.  Calculated  values  of  the  amplitude  (a)  and  phase  (b)  coherence 
lengths  as  functions  of  the  propagation  angle  and  source  coherence 
length  (pe)-  Fixed  parameters  are  X  -  10“®  m,  -  10"''®  L  - 

10"''  m,  and  =  10®  m. 

on  the  predicted  field  coherence,  Fig.  4  shows  calculated  values 
of  the  amplitude  and  phase  coherence  as  functions  of  the 
propagation  angle  and  source  coherence  (for  relatively  inco¬ 
herent  sources),  while  Fig.  5  depicts  the  same  quantities  as 
functions  of  the  strength-of-turbulence  (€„)  and  source  co¬ 
herence  for  nearly  coherent  sources).  Two  different  illus¬ 
trations  are  employed  because  the  vacuum  intensity  angular 
distribution  [Eq.  (44)]  becomes  strongly  peaked  in  the  axial 
direction  for  very  coherent  sources  and  hence  angular  plots 
of  the  predicted  coherence  lengths  become  meaningless  for 
relatively  coherent  sources.  It  is  evident  from  Fig.  4  that  the 
angular  distribution  of  the  amplitude  coherence  length  [Fig. 
4(a)]  is  relatively  unaffected  by  the  source  coherence  length 
(for  moderate  turbulence  and  the  range  of  coherence  lengths 
plotted)  although  the  phase  coherence  length  decreases  in 
proportion  to  the  inverse  square  of  the  source  coherence 
length.  Because  of  the  nearly  rectilinear  propagation  char¬ 
acteristics  of  coherent  sources  Fig.  5  shows  the  phase  and 
amplitude  coherence  lengths  for  nearly  coherent  sources  at 
propagation  angle  of  10"^  rad  and  0®,  respectively,  (recall  that 
^ir/4  as  0°)  as  a  function  of  the  refractive-index 
structure  constant  and  the  source  coherence  length.  The 
/ 

182  J.  Opt.  Soc.  Am.,  Vol.  68,  No.  2,  February  1978 


log,o  [c2(m-2/3)] 


Pe<m) 


FIG.  5.  Calculated  values  of  the  amplitude  (a)  and  phase  (b)  coherence 
lengths  as  functions  of  the  refractive-index  structure  constant  (C^)  and 
source  coherence  length  (pe)-  Fixed  parameters  are  X  =  10"®  m,  /.  =  10"'' 
m,  R'^  -  10®  m,  and  0'^  -  0®  [Fig.  5(a)]  and  ==  10"®  rad  [Fig.  5(b)]. 

decrease  of  the  predicted  phase  coherence  length  with  in¬ 
creasing  source  coherence  is  again  illustrated  in  Fig.  5(b),  and 
it  is  additionally  noted  that  the  strength-of-turbulence  affects 
this  functional  dependence  appreciably  only  for  the  highest 
refractive-index  structure  constants  (where  the  atmospheric 
coherence  length  becomes  comparable  to  the  source  coherence 
length).  Significantly,  the  amplitude  coherence  length  in¬ 
creases  rapidly  for  the  high  range  of  source  coherence  lengths 
plotted  in  Fig.  5(a),  whereas  the  turbulence  strength  influ¬ 
ences  this  functional  dependence  at  the  high  values  of  the 
refractive-index  structure  constant. 

A  final  illustration  of  the  functional  dependence  of  the  field 
coherence  predicted  by  the  above  analysis  is  provided  by  the 
Fig,  6  plots  of  the  amplitude  and  phase  coherence  lengths  as 
functions  of  the  source  size  (L)  and  refractive-index  structure 
constant.  Both  phase  and  amplitude  coherence  lengths  dis¬ 
play  a  similar  functional  form  with  respect  to  the  variables 
plotted  in  Fig.  6.  Clearly,  the  amplitude  coherence  lengths 
become  speckle  dominated  for  weak  turbulence  (Cl  <  10"®) 
as  evidenced  by  the  constant,  asymptotic  relationship  with 
respect  to  the  strength-of-turbulence  variable  (Cl).  Simi- 
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larly,  the  phase  coherence  length  is  constant  for  weak  turbu¬ 
lence.  However,  in  each  case,  atmospheric  effects  reduce  the 
respective  coherence  lengths  at  high  refractive-index  struc¬ 
ture-constant  values  with  the  onset  of  the  reduction  occurring 
earlier  for  the  smaller  source  dimensions.  The  data  plotted 
in  Fig.  6  display  trends  similar  to  the  predictions  of  Lee  et  al. 
for  an  incoherent  source  in  the  axial  direction  (Figs.  3-5  of  Ref. 

'7)- 

REMARKS  AND  CONCLUSIONS 

The  theoretical  foundation,  provided  above,  has  several 
significant  advantages  over  previous  investigations  of  tur¬ 
bulence  effects  on  radiation  propagation.  The  principal 
contributions  of  this  analysis  are  (i)  it  is  mathematically 
straightforward,  (ii)  it  provides  for  the  characterization  of 
radiation  derived  from  partially  coherent  sources  throughout 
the  full  range  of  propagation  directions,  (iii)  it  includes  the 
known  features  of  (a)  coherent  radiation  propagation  in  tur¬ 
bulence  and  (b)  partially  coherent  radiation  propagation  in 
vacuo,  and  (iv)  it  permits  a  physical  characterization  of  the 


FIG.  6.  Calculated  values  of  the  amplitude  (a)  and  phase  (b)  coherence 
lengths  as  functions  of  the  source  size  (L)  and  refractive-index  structure 
constant  (C^).  Fixed  parameters  are  X  =  1 0“®  m,  pe  =  10“'*  m,  ^  =  1 0^ 
m,  and  0+ =  5®. 
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radiation  coherence  propagation  in  terms  of  scattering  func¬ 
tions.  Although  the  five-thirds  law  structure  function  was 
approximated  by  a  quadratic  function  to  permit  explicit  cal¬ 
culations  of  coherence  effects,  calculated  plane-wave  results 
are  in  sufficient  agreement  with  experimental  data  to  permit 
confidence  in  the  general  trends  and  magnitudes  of  coherence 
effects  predicted  by  this  approximation.  Significantly,  all  of 
the  predicted  coherence  effects  derived  from  this  approxi¬ 
mation  are  consistent  with  physical  intuition  and  approach 
proper  limiting  results  over  an  appropriate  range  of  propa¬ 
gation  angle,  source  size,  source  coherence  length,  and  tur¬ 
bulence  strength  variables.  Although  the  definition  of  an 
effective  far-field  range  for  Gaussian  partially  coherent 
sources  resulted  from  considerations  unrelated  to  turbulence 
effects,  this  definition  is  significant  because  it  extends  the 
concept  of  the  far  field  to  partially  coherent  sources.  While 
this  extension  is  limited  to  partially  coherent  sources  with 
Gaussian  MCF’s,  it  is  likely  that  similar  results  obtain  for 
other  source  MCF’s  because  of  the  Fourier  transform  rela¬ 
tionship  [Eq.  (20)]  between  the  field-point  MCF  and  the 
generalized  source  coherence. 

The  above  analysis  of  the  atmospheric  propagation  of 
partially  coherent  radiation  should  find  application  in  areas 
where  laser-illuminated  objects  viewed  through  the  atmo¬ 
sphere  are  of  interest  (e.g.,  laser  radars).  Analyses  have  been 
developed  (e.g..  Refs.  25  and  26)  for  predicting  the  mutual 
coherence  and  intensity  distribution  of  laser-illuminated 
rough  surfaces.  Because  the  foregoing  development  permits 
analysis  of  atmospheric  modifications  of  vacuum  distributions 
[via  the  convolution  Eqs.  (24)  and  (32)],  theoretical  investi¬ 
gations  of  turbulence  effects  on  realistic  laser-illuminated 
rough  surfaces  should  be  possible. 
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APPENDIX  A 

Phase  factor  evaluation 
The  phase  factor  defined  in  Eq.  (6)  can  be  written 

L(r+,  r~,  P+,  P-)  =  |r+  -F  -  P+  -  %P-|  -  |r+  -  %r“ 

-  P+  +  VsP”!  =  [|r+  -  P+|2  -f  (r+  -  P+) .  (r“  -  p-) 

+  ii4|r“  “  P-|2]i/2  ^  [|r+  _  p+|2  ^  __  p+) .  (y-  _  p-) 

+  y4|r-  -  p-|2]i/2  =  |r+  _  p-f-m+y2[(r+  _  p+)..  (r-  p-) 

+  y4|r“-P“|2]|r+  -P+|-2-F0(|r+  -P+|“2)  i 
-  V2[-(r+  -  P+) .  (r-  -  P-)  -F  lAIr-  -  P-|2]|r+  -  P+|2 
-0(lr+--P+|-2j^  (r--p-).t£)(r+  p+)^  (ai) 

where 

w(r+  P+)  H  (r+ -  P+)/|r+ -  P+l 

=  (r+  -  +  (r+)2(i?+)-2 

-2(r+.s+)(ft+)-2]-i/2  (A2) 

and  the  definitions  of  Eqs.  (7)-(10)  have  been  employed.  The 
vector  u;(r+,  P+)  is  illustrated  in  Fig.  1.  Because  is  much 
greater  than  a  source  dimension  (by  assumption), 
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(A3) 


where 


M;{r+  P+)  s  (r+/R+-s+). 
Using  Eq.  (A3)  and  the  fact  that 

P-  =  R+s-  +  R-S+, 


(A4) 


together  with  the  knowledge  that  s"*"  and  s  are  orthogonal 
vectors  (s+  is  a  unit  vector),  Eq.  (Al)  can  he  written 

L(r+,  r-,  P+,  P-)  =  [r-  •  r+/R+  -r--s+-r+-s- 

+  R-(l-^+T+/fl+)].  (A5) 

However,  again  noting  the  small  value  of  (r+/f?+)  over  the 
range  of  integration  considered,  Eq.  (A5)  becomes 

L(r+,  r-,  P+,  P-)  s  [r-  •  r+/fl+  -  r"  •  s+  -  r+  •  +  ft"] 

(A6) 


to  an  excellent  approximation. 


APPENDIX  B 

Uniform  coherent  source  MCF  and  spectral  density 
evaluation 

A  uniform  coherent  source  is  characterized  by  a  source 
MCF  of  the  form 

r,(r+  r-)=  |Uol2.  (Bl) 

Therefore,  the  generalized  two  field-point  source  coherence 
[Eq.  (15)]  is  given  by 

G  =  (2ir)2(/J+/fe)2A-i|Uo|25(r-  -  «+s*-),  (B2) 

and  hence  the  predicted  MCF  is  [Eq.  (17)], 
rturb(p+^  P-)  =  I  Uo|2  C0s2fl+ 

X  J  J'”  d2i.-e-i(fer-H+5(r  -  R+i-)F(t-.  P") 

=  I  UoP  cos2e+e'*«^*‘-^*F(/?+s -,P-) 

=  |17o|2cos20+F(P-,P-),  (B3) 

assuming  equidistant  radial  field-point  separations  (i.e.,  R~ 
=  0).  Equation  (38)  can  be  easily  integrated  using  the  con¬ 
ditions  of  equal  source-point  and  field-point  separations  im¬ 
posed  by  Eq.  (B3);  the  result  is 

F(P”,  P-)  =  exp[-3(p/po)2],  (B4) 

where 

P=1F(P-)|. 

Equations  (B3)  and  (B4)  express  the  plane-wave  MCF  result 
using  the  quadratic  approximation  [Eq.  (38)].  Note  that  the 
coherence  length  predicted  by  this  approximation  is 

p(e->)|plane.wave  =  (l/V3)po  =  0.58po,  (B5) 

which  is  in  good  agreement  with  the  value  of  0.56  po  predicted 
by  the  five-thirds  law. 

Artem’ev  and  Gurvich^i  have  performed  optical  measure¬ 
ments  of  the  spectral  density  of  field  fluctuations  using  the 
Fourier-transform  properties  of  a  simple  lens.  Accordingly, 
their  dimensionless  spectral  density  results  are  defined  in  the 
following  terms: 

V(x)=c/(x)/m  (B6) 


J(y)  =  r  “  2Rot)  cos[{2kRoy/Fo)t]  dt,  (B7) 

Jo 

X  =  2kRoyFonDi(2Ro)]-‘^^^  (B8) 

Fo  is  lens  focal  length,  2Ro  is  lens  diameter,  and  Di  is  the 

complex  phase  structure  function  for  the  assumed  plane  wave. 
Fourier  transforming  the  plane-wave  MCF  predicted  by  Eqs, 
(B3)  and  (B4)  according  to  the  prescription  of  Eq.  (B7)  yields 
the  following  result: 

y(xO  =  exp  [-(l/48)x'^p^o^Di(2Fo)],  (B9) 

where  the  dimensionless  variable  x  has  been  redefined 

x'  =  2kRoyFo^[Dii2Ro)]^^^^  (BIO) 

in  accord  with  the  quadratic  assumption.  Because  the 

plane-wave  MCF  result  in  the  quadratic  approximation 
implies  a  structure  function  defined  by 

V2Z)i(p)  =  3(/)/po)2,  .  (BID 

Eq.  (B9)  reduces  to 

V(x)  =  exp(-xV2).  (B12) 

Equation  (B12)  is  plotted  in  Fig.  2  together  with  experimental 
data  of  Artem’ev  and  Gurvich.21 

APPENDIX  C 

Two-source  spherical  wave  MCF  evaluation 

To  perform  the  projection  operations  indicated  in  Eq.  (38), 
orthogonal  unit  vectors  di  and  02  are  constructed  which  define 
a  plane  transverse  to  the  P"*”  direction.  These  vectors  are 

.di=[s+X/e][l-(s+.^)2]-l/2 

=  —I  cos<t>  ^  j  sin0,  (Cl) 

02  =  Oi  X  S”*" 

=  — f  cos^"^  sin0  +  j  cos8'^  cos0  +  sin^'^^.  (C2) 

Therefore, 

^(r")  =  (di  •  r")di  +  (02  •  r”)d2 

=  —di(x~  cos(^  +  y”  sinij)) 

+  d2{—x  ~  sin<^  +  y~  cos(/>)  cos0  (C3) 

The  vector  difference  vector  can  be  written 

§~  =  §+  X  dfl,  (C4) 

where 

s'^  =  i  sin4>  sin^'*'  ~  j  sin^'*'  cos</>  +  ^  cos^"*",  (C5) 

and 

dQ  ^  6  COS0/  +  8  sixKpj  -h  (C6) 

is  the  axial  vector  corresponding  to  the  antisymmetric  rotation 
tensor  that  carries  s'  into  s  (in  the  small-angle  approximation). 
Performing  the  vector  multiplication  of  Eq.  (C4)  using  the 
expressions  (C5)  and  (C6)  and  performing  the  projection  op¬ 
eration  indicated  in  Eq.  (C3)  yields 

F(P“)  =  F‘^(a>  sin^di  +  602),  (C7) 
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where 


P--/?+s-  +  /?“s+.  (C8) 

Because  the  problem  is  azimuthally  symmetric,  the  plane 
included  by  the  y  axis  is  chosen  (i.e.,  0  =  tt)  so  that 

\VP{F-)  +  (1  -  tyP{r-)\^  =  t^(R+o)  sin0+  -  a:-)2  +  (i?+6 
+  y~  cos0'^)2]  +  2t[x~(R'^(jo  sinO'^  ~  jc”)  —  y“  cosO'^iR'^d 
+  y'^  cos^"*")]  +  [(jc“)2  +  (y“  cos^'^)^].  (C9) 

Employing  the  definitions  of  Eqs.  (48)  and  (53)  yields 

Fir-,  P",  R^)  =  exp  \-(R/pon{b  + 

+  (c  —  5)2  +  c6  “  6d]j  (CIO) 

upon  integrating  Eqs.  (38)  using  Eq.  (C9). 

Utilizing  the  definition  of  the  pseudo-scattering  function 
[Eq.  (22)]  together  with  the  above  result  for  the  two-source 
spherical-wave  MCF  yields 

F(hs-,R+)  =  (fe«+)2sec0+  exp[-^(— )V  +  82)j 

X  J  |-  (^t(‘  *  ('  -011 

X  exp[--  i{kR‘^)(^xb  +  sec^+^^c)] 

=  7r(A:po)^  sec^+  exp  -  -  (~)  ] 

X  exp  (-  i  ikpo)^  +  ^2  sec20+))  ] 

X  exp  [i  -  ^yd  secfl+)]  (Cll) 

after  changing  variables  and  integrating. 
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